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Abstract 

In hadron induced air showers of highest energies (E>10 18 eV), the 
lateral distribution functions of electrons and muons are a superposition 
of many separate electromagnetic sub-showers, initiated by meson decay 
at different altitudes and energies. The lateral distribution function is 
the primary tool for reconstructing the energy of the primary in a UHE 
cosmic ray shower, so understanding it in detail is a prerequisite for having 
confidence in the energy determination. We analyze in this paper the 
dominant contributions to the ground level lateral distribution functions, 
as a function of the altitude and energy at which the sub-showers are 
initiated. Far from the core, the dominant contribution to the density 
of electrons comes from sub-showers initiated at low altitudes and low 
energies (E<100 GeV). The dominant sub-showers are initiated at large 
radial distance from the core and at a large angle with respect to the 
main shower axis. This demonstrates the need for careful treatment of 
low energy hadron physics modeling even for ultrahigh energy primaries. 



1 Introduction 

There are basically two methods to measure air showers induced by cosmic rays 
with an energy greater than 10 18 eV. One is to measure the fluorescence light 
caused by charged particles passing through the atmosphere; the other is to 
sample the lateral distribution of charged particles on the ground. Theoretical 
predictions for both the longitudinal and lateral distributions of secondaries are 
produced by Monte Carlo simulation of the showers. Comparing theory and 
observation permits properties of the primary cosmic ray, such as its energy 
and composition, to be inferred. 
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In this paper we study lateral distribution functions of electrons and muons. 
In particular we want to determine which altitudes and which regimes of energy 
of the hadronic part of an air shower are of primary importance in producing 
various components of the ultimate ground level lateral distribution. This is a 
step toward the larger goal of elucidating the sources of systematic uncertainties 
in theoretical predictions for shower properties. 

2 Air showers 

Atmospheric showers can be categorized into three elements, the hadronic part, 
the electromagnetic part and the muonic part. Basically, the hadronic part feeds 
the other two components, photons being mainly produced by ir /rj decay and 
muons by the decay of charged pions and kaons. Once a photon is produced it 
initiates an electromagnetic sub-shower which evolves essentially independently 
of the rest of the shower. (With some low probability, a photon interacts hadron- 
ically in a photo-nuclear reaction, re-initiating a hadronic sub-cascade; this is 
neglected here.) Muons, once created, interact very little. For them the most 
important process is energy loss via ionization, and occasionally pair production 
or bremsstrahlung in the nuclear field of air molecules. 

Let us call the photons and muons which are directly produced by hadronic 
decays, initiators, because they initiate a sub-shower of their kind. The ques- 
tions we address here are: where are these initiators produced, what are their 
energies, and what are the energies of the hadronic reaction in which they were 
themselves created. This is crucial to understanding lateral distribution func- 
tions and especially what governs the particle densities at large distances from 
the shower core. 

For electromagnetic sub-showers there are two classes of initiators: high 
energy photons produced at high altitudes and low energy photons produced 
at low altitudes because high altitude initiators must have sufficient energy for 
their secondaries to reach ground level before being absorbed and sub-showers 
initiated at low energy must be created at low altitude, to avoid being absorbed 
before reaching the ground. However the the relative importance of these classes 
of initiators, as a function of the radial distance of the observed particle from 
the sub-shower, and the dominant sources of lateral spreading and arrival-time 
delay, has not been systematically studied until now. 

For muons the situation is simpler than for electromagnetic showers. High 
energy charged pions rarely decay due to their huge Lorentz factor, so only 
low energy components of the hadronic shower contribute significantly to muon 
production. Moreover Lorentz beaming means that lower energy particles have 
typically larger angles with respect to the shower axis, and therefore contribute 
more to the muon density at larger distance from the core, than do higher energy 
particles. For a given initiator energy and angle, simple geometry implies that a 
muon can propagate to larger distances from the core if it is produced at higher 
altitudes. Muons, in contrast to electrons and photons propagate without much 
absorption to the ground level. 
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The air showers analyzed in this paper are simulated with a traditional 
Monte Carlo method, whose code will be released soon. The basic features of 
this model are the high energy hadronic model QGSJET [Q, the low energy 
hadronic model GHEISHA ||] and the electromagnetic shower model EGS4 
[||. The model has also some new simulation techniques which speed up the 
computation of air showers considerably, but these are irrelevant for the present 
analysis. Therefore, the pure physics content (hadronic and electromagnetic 
modeling) is identical to the most commonly used option in the CORSIKA 
package [Q. We have directly confirmed the agreement of the results presented 
here with CORSIKA simulations. 

3 Initiators 

In the following we analyze the contributions to the lateral distribution functions 
of 5 x 10 19 eV proton-induced vertical showers at a ground level altitude of 900m, 
averaging over 50 showers. The kinetic energy cutoffs for registering hadrons, 
muons, electrons/positrons and photons have been chosen to 0.05 GeV, 0.05 
GeV, 1 MeV and 1 MeV, respectively. 

3.1 Ground electrons 

For electrons at ground level, the results are best summarized in Fig. [l], which 
shows the relative contribution at three distances from the core as a function 
of the electromagnetic initiator's energy and production depth. This is just the 
probability that an electromagnetic sub-shower is initiated at a given altitude 
and energy. The altitude is plotted in units of vertical atmospheric depth, 
defined as X(h) = f. p(h')dh', h being the altitude and p(h) the density of 
air, according to a standard USGS atmosphere. The results are obtained as 
follows: for each electron arriving at a certain distance from the shower core 
we trace back the history of this particle, generation by generation, and record 
the energy and depth of production of the initiator (the last particle of the 
sub-shower which was not a hadron: in most cases, a photon from ir°/r] decay). 

One recognizes two distinct regions of initiator energy and depth of pro- 
duction which play the dominant role in producing electrons, depending on the 
distance from the core. Closer to the core, the initiators come predominantly 
from high altitudes and have high energies, while far from the core the initiators 
are seen to be primarily produced at low altitude with low energy. 

We can analyze this interesting behavior in greater detail by looking at the 
probability distributions as a function of energy and depth separately. Figure 
g shows the distribution of energies of electron initiators whose electromagnetic 
sub-showers contribute to the electron density at 4 different lateral distances. 
The area under each curve is normalized to 1 to get probability distributions. 
Close to the shower core the most important contributions come from high 
energies. At larger distances the contribution of low energy initiators becomes 
more and more important. Fig. || gives the probability distribution as a function 
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of the altitude of the initiators. At large distances from the core the dominant 
contribution comes from relatively low altitudes, i.e. large depth. This confirms 
that the power of lateral spreading for each high energy electromagnetic cascade 
is limited. In fact, the greatest contribution at large lateral distance comes from 
low energy sub-showers which are emitted at a large angle with respect to the 
main shower axis, as seen in Fig. Q and at a large distance with respect to the 
shower axis (Fig. ||). It is useful to read Figs. |^ and || together, to see at which 
altitudes the initiators are produced, and then to read off the mean distance 
from the shower axis. One realizes that high altitude initiators are produced 
close to the axis, and low altitude ones close to the point where the electrons 
reach the ground. This means that the low energy hadronic part of the shower 
is producing a significant portion of the lateral expansion. 

Fig. H shows a very irregular structure below 300<?/cm 2 . This is a statistical 
fluctuation caused by some events having a very energetic photon, which is 
consequently a sub-shower initiator for many electrons and photons at ground 
level. 

We note that the results on electron initiators presented here suggest why 
analytic treatments of the lateral distribution function of the air showers of 
hadronic primaries using a formula like the NKG M, |(| are unsatisfactory, even 
though the formula may perfectly describe an electromagnetic sub-shower. The 
implementation of such analytic approaches generally assumes that each sub- 
shower initiator has a small angle with respect to the shower axis. As we have 
seen, this is not correct when describing electrons and positrons at large dis- 
tances from the shower core. Only an implementation of NKG with full consid- 
eration of the emission angle of each sub-shower initiator, without small angle 
approximation, and taking the correct intersection of the inclined sub-shower 
with the ground level, can be expected to give reasonable results. 

3.2 Ground photons 

Photons show similar behavior to electrons and positrons, since they are created 
in the same electromagnetic sub-shower. However, Figs. | and | reveal that the 
effect discussed above for electrons is much less pronounced. Photon initiators 
are produced higher in the atmosphere than electron initiators, even though 
they are part of the same electromagnetic sub-showers. This is probably due to 
the fact that photons do not lose energy by ionization loss, therefore low energy 
photons have a higher chance to propagate a large distance than electrons do. 

3.3 Ground muons 

The corresponding results for muons can be seen in Figs. || and||. Their energy 
spectrum (Fig. ||) at the point of creation gets peaked at a lower energy as 
one selects muons arriving at larger distance from the shower axis, because a 
lower energy makes larger emission angles more probable. As for the depth of 
creation, one can see in Fig. || that at high altitudes muon initiators contribute 
more at large distances from the core, but never to the extent of electromagnetic 
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initiators. This is not surprising since the hadronic shower has to evolve down 
in energy in order for pions or kaons to decay. 

A comparison of Figs. || and [To] shows that muons observed at large radial 
distance are preferentially but not exclusively produced close the the shower 
axis, propagating to larger distances due to their large angle of production with 
respect to the shower axis. 

The small peak in Fig. |9| at r=2000m and very low altitudes may be a 
fluctuation or the onset of a low altitude contribution as observed for electron 
densities; however the phenomenon is much less pronounced and we do not 
investigate it further here. 

4 Hadronic reaction energy 

It is also interesting to investigate the energy of the hadronic reactions which are 
responsible for creating the sub-shower initiators. Figure [ll] shows the energy 
of the projectile in the lab frame, of the hadron-air collision which creates the 
electron- or muon-initiator. It reveals that high energy hadronic physics governs 
the production of the electromagnetic component of the ground level shower 
close to the core, but for larger distances from the core the last hadronic reaction 
is in the low energy (E<100 GeV) regime. 

Figure [l^ shows that at large distances from the core, low energy reactions 
become increasingly important for determining particle densities. The top figure 
shows the total lateral distribution function of electrons plotted in a full line, 
while the part contributed by initiators which are created in an hadronic reaction 
with an energy > 100 GeV is shown with a dashed line. The corresponding 
curves for muons are drawn as dotted (muon lateral distribution function) and 
dashed-dotted lines (contribution by high energy reactions with E\ a h > 100 
GeV only). The bottom figure shows the ratios of the corresponding curves, for 
electrons and for muons. At 600m 20% of the electron density comes from a 
low energy hadronic reaction; at 1km this fraction is 40%. The situation is even 
more pronounced in the case of muons, due to the comparatively long lifetime 
of charged mesons. Here the last hadronic reaction is in most cases a low energy 
reaction, even at small distances from the core. 

These results shed new light on the importance of low energy hadronic mod- 
eling in air showers, an aspect which is often neglected in the literature. The 
impact of low energy shower modeling can be directly seen in Fig. [□]. In the sim- 
ulation program used for this study, we employed the low energy hadronic model 
GHEISHA to describe interactions of projectiles whose lab energy is less than 
90 GeV, while at higher energies QGSJET is used. This follows the CORSIKA 
procedure, only differing by choosing the energy of transition between low and 
high energy hadronic models at 90 GeV rather than 80 GeV. Fig. [ll] shows that 
the projectile energy distributions show a break at log 10 (-E/GeV) « 2 , espe- 
cially for muons (lower figure) at large distance from the core. The sharp break 
is an artifact of a discontinuity in the modeling of the hadronic interactions, 
between QGSJET and the low energy hadronic model GHEISHA. Evidently, at 
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the transition energy E\ a h — 90 GeV, GHEISHA produces many more muons 
at large distances than QGSJET; this arises because the two models have dif- 
ferent differential cross sections and particle production levels. (For examples 
of differences between these models in this energy region, see With correct 
modeling, the high and low energy models must match at the transition energy, 
in principle in every detail. What we see here is that at large radius, these 
details affect the predicted shower properties. We have checked, however, that 
the qualitative conclusion that low energy, low altitude initiators dominate for 
large radius electrons, is not an artifact of the non-smooth transition between 
low and high energy hadronic models. 

It has already been pointed out that for much lower energy primaries the 
low energy hadronic model is important for muons measured at ground level. 
In Ref. |g, ^| it is shown that for 10 15 eV hadronic showers, the energy of the 
last hadronic reaction before producing a muon (E^ > 1 GeV) is in 50% of all 
cases a low energy reaction with E\ a ^ < 100 GeV. Refs. Jl(], Q study the low 
energy hadronic model dependence by analyzing CORSIKA |Q| simulations of 
10 15 eV hadron induced showers with two different models at energies below 



80 GeV: GHEISHA|§ and UrQMD (12[ |13|. The two models show differences 



in the energy spectrum of produced muons, but the overall number of muons 
stays approximatively the same in both cases |Tl[ . This fact becomes important 
when measuring muons with different energy thresholds as shown in Ref. ||l0| . 
While possibly not surprising, the analysis presented here shows that the low 
energy hadronic model is also important for determining the lateral distribution 
function at large radius for ultrahigh energy cosmic ray showers. 



5 Mean arrival time of ground particles 

We have seen that in the case of lateral distribution functions of electrons, there 
are two distinct regions in altitude and energy which contribute differently as a 
function of distance from the shower axis. This gives rise to the question where 
the time delay of particles with respect to the shower front comes from. In a 
simplified model one can consider three contributions: The time delay is caused 
by 

• hadronic showering, 

• electromagnetic showering, and 

• geometry. 

Fig. |l^ shows illustrates the different parts. The hadronic time delay is defined 
as Aihad = tAB — AB/c with AB being the distance from A to B and tAB 
being the actual time spent in the hadronic cascade prior to production of the 
initiator, i.e. going from point A (first interaction) to point B (production of the 
initiator, e.g. a photon or a muon). Ai cm = tsc — BC/c is the corresponding 
time delay of the electromagnetic (or muonic) sub-shower. The total time delay 
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is defined as A< to t = tAB +tBC — AC/c. And finally the geometric time delay is 
Ai gco = (AB + ~BC - AC) / c = At to t - Ai had - Ai em . A detailed Monte Carlo 
model permits keeping track of each propagating particle and hence computing 
the values Ai to t, Ai had , Ai Gm . 

Fig. [l^ shows the results, which are best understood by referencing also 
the figures of the previous sections giving production height and radius of the 
initiators. Ground electrons close to the core spend most of their time in the 
electromagnetic sub-shower, since for them the main contribution comes from 
high altitude high energy initiators; the hadronic part of their history is at high 
energy for which multiple collisions do not increase the path length significantly, 
so their hadronic time delay is small. Further away from the core the hadronic 
part of the time delay increases, becoming as important as the electromagnetic 
one at 1500m. The increasing hadronic time delay is correlated with an increas- 
ing time delay caused by geometry, since both the radius of production of the 
initiators and their average emission angle increases with radius of the observed 
electron. Thus at very large distances, geometry and hadronic time delay are the 
dominant contributions, and the electromagnetic spreading is roughly constant 
as we have seen earlier. 

As for ground photons, the main time delay is caused by the electromagnetic 
sub-shower. The difference from electrons becomes clear when one remembers 
that photon initiators are higher in the atmosphere than electron initiators (see 
Figs. | and @). 

The time delay of muons is mainly caused by geometry and, though less 
important, the hadronic part of the shower. 

6 Conclusions 

We have analyzed in detail how electromagnetic and muonic sub-showers con- 
tribute to the overall lateral distribution functions of an ultrahigh energy cos- 
mic ray air shower. At all distances from the core, ground level muons come 
mostly from low energy muons emitted at intermediate-high altitudes in the 
atmosphere. However, at large distances from the core (>1000m), electrons 
come predominantly from low energy n°/r/ 's produced at low altitudes and 
large distance from the core, and with a large angle with respect to the shower 
axis. Moreover the collisions producing the initiators are themselves low en- 
ergy hadron-nucleus interactions which in standard air shower simulations are 
treated by a different model than the high energy interactions. 

We conclude that in order to achieve detailed agreement between data and 
models for the lateral distribution functions and arrival-delay time distributions, 
necessary in order to optimally reconstruct properties of the primary, it is nec- 
essary to take greater care for modeling the low energy hadronic interactions in 
the showers. 
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Figure 1: The relative contribution to electron densities at three distances from 
the core r<100, r=1000, and r=2000m, as a function of energy and production 
depth of the initiator of the electromagnetic sub-shower. Two distinct regions 
are visible, one at low energy and low altitude, and another at high altitude and 
high energy. 
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Figure 2: The top figure shows the probability distribution of initiators for 
ground level electrons as a function of the initiator's energy for different radial 
distances from the core. The bottom figure shows the integrated probability 
P(E' < E) = f Q j^ydE' which corresponds to the fraction of initiators with 
energy less than E. At r=2000m 65% of the electron density comes from low 
energy gammas with E' < 1 GeV. 
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Figure 3: The depth of electron initiators for 4 different distances from the 
shower core. The greater is the distance from the core, the larger the contribu- 
tion from low altitudes. 
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Figure 4: The emission angle of electron initiators as angle with respect to 
the main shower axis. The first bin has been duplicated in negative values for 
better visibility. Close to the core, the initiators are predominantly emitted 
collinear with the shower axis, while at larger distances greater emission angles 
are dominant. 
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Figure 5: The mean distance of electron initiators to the shower axis. 
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Figure 6: Energy distribution of photon initiators. 
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Figure 7: The distribution of vertical depth of photon initiators. 
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Figure 8: This plot shows the energy distribution of muon initiators, which 
the energy of muons at the point of creation, leading to ground level muons 
different lateral distances. 
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Figure 9: The depth of muons at the point of creation. 
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Figure 10: The mean distance of muon initiators to the shower axis. 



15 




energy logio(E/GeV) 




energy logio(E/GeV) 

Figure 11: Energy of the last hadronic reaction before creating a electromagnetic 
(top) and a muonic (bottom) sub-shower. 
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Figure 12: The importance of high energy reactions to the density of electrons 
and muons as a function of distance from the shower core: the full line in the 
top figure shows the total density of electrons as a function of distance from the 
shower core. The dashed line represents the density of electrons whose initiators 
are created in a high energy reaction (-Ei a b > 100 GeV) . The corresponding 
curves for muons are plotted as dotted and dashed-dotted lines. The bottom 
figure shows the ratios for electrons (full line) and muons (dashed line). 
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Figure 13: The contributions to the time delay of particles: hadronic time delay 
from A to B, electromagnetic/muonic from B to C. The difference of the sum 
(dotted line) to the total time delay (dashed line) represents the geometric delay 
as explained in the text. 
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Figure 14: Contributions to the mean time delay of electrons (top), photons 
(middle) and muons (bottom) as a function of the distance from the shower 
axis. 
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